Draft version September 28, 2009 

Preprint typeset using I^T^X style cmulateapj v. 08/22/09 



KINETIC SIMULATIONS OF TURBULENT MAGNETIC-FIELD GROWTH BY STREAMING COSMIC RAYS 
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ABSTRACT 

Efficient acceleration of cosmic rays (via the mechanism of diffusive shock acceleration) requires 
turbulent, amplified magnetic fields in the shock's upstream region. We present results of multi- 
dimensional particle-in-cell simulations aimed at observing the magnetic field amplification that is 
expected to arise from the cosmic-ray current ahead of the shock, and the impact on the properties 
of the upstream interstellar medium. We find that the initial structure and peak strength of the 
amplified field is somewhat sensitive to the choice of parameters, but that the field growth saturates 
in a similar manner in all cases: the back-reaction on the cosmic rays leads to modification of their 
rest-frame distribution and also a net transfer of momentum to the interstellar medium, substantially 
weakening their relative drift while also implying the development of a modified shock. The upstream 
medium becomes turbulent, with significant spatial fluctuations in density and velocity, the latter in 
particular leading to moderate upstream heating; such fluctuations will also have a strong influence 
on the shock structure. 

Subject headings: acceleration of particles, cosmic rays, methods: numerical, shock waves, supernova 
remnants, turbulence 



1. INTRODUCTION 

The forward shocks of young shell-type supernova rem- 
nants may be efficient acceleration sites for Galactic cos- 
mic rays. The theory of diffusive shock acceleration 
(DSA) provides a promising mechanism by which a small 
fraction of particles can attain nonthermal energies in 
such an environment, provide d they are c onfine d to the 
shock vicinity; the review by iRevnoldsl (|2008h and its 
references have treated the subject in detail. This con- 
finement may be supplied to some degree by the inter- 
stellar magnetic field, but the estimated upper limit in 
energy falls short of what is observed in cosmic rays of 
Galactic origin. The implied need for a much stronger, 
turbulent upstream field, together with observational ev- 
idence for strong magnetic fields immediately behind the 
forward shock, has attracted attention to the question of 
upstream magnetic-field amplification. 

It is well known that ion beams can generate mag- 
netic t urbulence via resonant a n d non - reson ant interac- 
tions (jWinske and Lerovl [l98l . [Bel (|2004h (hereafter 
B04) suggested that non-resonant magnetic field ampli- 
fication may also be caused by the cosmic-ray current 
expected in the cosmic-ray precursor of a quasi-parallel 
shock. Although substantial field amplification has oc- 
curred in magnetohydrodynamical (MHD) simulations 
that assume a constan t cosmic-ray current (B04; also 
iZirakashvili et all [2008), the approximations of MHD 
may not be appropriate for modeling the non-linear evo- 
lution and eventual saturation of this instability; we must 
use kinetic methods to simulate these later stages with 
accuracy, including any back-reaction on the cosmic rays. 

1 Now at: Institut fur Physik und Astronomie, Universitat Pots- 
dam, 14476 Potsdam-Golm, Germany; 
and DESY, 15738 Zeuthen, Germany 



Our earlier work (Niemiec "eTaLl [20081 hereafter N08) 
used 2-D and 3-D particle-in-cell (PIC) simulations to 
model the growth and saturation of a current-driven in- 
stability. We relaxed the condition assumed in the cal- 
culations of B04 that the rate of unstable growth must 
be much less than the ion gyrofrequency, and we found 
that an oblique filame ntary mode dominated the initia l 
magnetic field growth. Riquclm e and Spitkovskvl (|2009D 
confirmed the findings of N08, and in particular inves- 
tigated the parameters for the transition between the 
filamentation mode seen in N08 and the parallel-wave 
mode seen in the MHD simulations, and verified that the 
parallel mode appears only in the regime with u> -C 
However, despite the difference in initial unstable modes, 
the non-linear cha racteristics of the system obser ved in 
both works (N08: [Riquelme and Spitkovskvl [20091 ) were 
similar. In particular, both investigations observed only 
modest amplification of the magnetic field, and satura- 
ti on by the same me c hanis m. The analytical calculations 
of ILuo and Melrose] (|2009f ). with a kinetic treatment of 
the current-driven instability, predict a saturation level 
consistent with the kinetic simulations. 

In this paper, we report more recent work, in which 
we simulate the current-driven instability and its satura- 
tion in a parameter regime in which the non-relativistic 
ion gyrofrequency far exceeds the (complex) frequency of 
the current-driven instability. Such parameters are par- 
ticularly relevant to the environment of non-relativistic 
strong shocks whose cosmic-ray precursors are sparse 
in comparison with the magnetized ambient interstel- 
lar medium (ISM); an understanding of the limits of 
field amplification in such an environment is therefore 
essential to the formulation of a self-consistent theory of 
cosmic-ray acceeleration in the Galaxy. We find that 
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the initial evolution of the magnetic field in this pa- 
rameter regime is in agreement with the predictions of 
B04 and cons istent with the transition parame ter, w/fij, 
identified bv iRiquelme and Spitkovskvl (|2009l ). We ob- 
serve that the field growth is limited to modest val- 
ues by the same saturation me chanism as before (N08; 
IRiquelme and Spitkovsk v 2009) . In addition to the mag- 
netic field, we examine the evolution of the upstream 
medium-precursor system in terms of bulk and local 
properties that are relevant for modeling the shock ac- 
celeration process. Notable findings include the rise of 
strong turbulent motion relative to the bulk upstream 
rest frame, and of motional electric fields associated with 
the magnetic fields that permeate the turbulent flow. 
The temperature of the upstream medium, measured in 
the local rest frame as determined with high spatial res- 
olution, attains significant values during the evolution- 
ary phase characterized by the turbulent motion, such 
that the ratio of magnetic to thermal energy density 
may in fact decrease as a result of this precursor-ISM 
interaction. Such modifications to the upstream medium 
may lead to substantial effects on the propagation of 
the subshock, and on the properties of the downstream 
medium: both the pre-existing upstream turbulence after 
advection and shock compression, and further magnetic- 
field generation via turbulent dynamo processes, may be 
quite sensitive to turbulence in the upstream medium 
(jGiacalone and Jokipu1l2007Hlnoue et al.ll2009[ ). 

We begin by describing in §2. li the simple physical sce- 
nario that captures all the ingredients necessary for the 
streaming instability we want to simulate, with i J2.2l sum- 
marizing our computational representation of that envi- 
ronment, and the specific parameterizations being out- 
lined in §2.31 We present our simulation results in $3] 
along with discussion of their significance, and finally 
summarize in f|4] 

2. SIMULATION SETUP 
2.1. Model 

To study the growth and saturation of magnetic turbu- 
lence in response to a cosmic-ray current, we simulate the 
evolution of a system representative of the upstream en- 
vironment of a non-relativistic strong parallel shock. We 
assume that the interstellar medium is a homogeneous 
collisionless plasma composed entirely of ions and elec- 
trons (number density iVj and N e , respectively) of oppo- 
site but equal charge ±e, permeated with a uniform mag- 
netic field £?|| . A population of cosmic ray ions (density 
Ncr — N e —Ni) drifts at the shock velocity v s h along the 
magnetic field relative to the plasma and carries a cur- 
rent density jcR = eNcRV s h \ the slight electron excess in 
the plasma maintains charge neutrality, and the electrons 
drift relative to the ions with Vd — v s hNcR/N e to provide 
a return current density j ret = —eN e Vd = —jcR- The 
ISM ion and electron populations are each distributed 
according to Maxwell-Boltzmann statistics in their re- 
spective rest frames with the same temperature, given 
by the electron thermal velocity v e _ t h- 

Two assumptions in our setup, that only electrons pro- 
vide the return current and that the streaming cosmic 
rays consist only of ions, arise from the large mass ra- 
tio rrii/m e between physical ions and electrons. The 
electrons will respond to any local charge or current im- 



balance on a much shorter timescale than the ions will. 
Moreover, at high energies, ions outnumber electrons (on 
account of their greater rigidity) by a factor that depends 
on the spectral index of their momentum distribution but 
is generally large enough that cosmic-ray electrons play 
a negligible role in the precursor physics. Cosmic rays 
upstream undergo pitch-angle scattering off magnetic in- 
homogeneities and approach a quasi-isotropic distribu- 
tion, forming a steady cosmic-ray precursor in the shock 
rest frame. Ions with lower energies will be advected 
back toward the shock sooner than ions with higher en- 
ergies, so at some distance upstream one expects to find 
a population of high-energy cosmic-ray ions drifting with 
the shock, up to some cutoff energy and/or distance 
(|Reville et al.ll2009h associated with particles that escape 
the system. 

For our simulation, we simplify the cosmic-ray distri- 
bution function. We assume that the isotropic, scattered 
cosmic rays of the precursor are energetically much more 
significant than any escaping cosmic rays that happen to 
be passing through the region of our simulation, so we 
neglect the latter altogether. Our box size is assumed to 
be small in comparison both to the scale length of the 
cosmic-ray precursor, such that we may neglect spatial 
variations in their distribution, and to the Larmor radii of 
any cosmic rays present, such that the details of their en- 
ergy distribution are similarly negligible, and we replace 
the already-narrow range of energies with a single energy 
(in the CR population rest frame), with Lorentz factor 
jcr- The quasi- linear calculations in B04 do not depend 
explicitly on cosmic-ray energy, and tests with various 
values of jcr have not produced significantly different 
outcomes, so it is unlikely that our choice of energy dis- 
tribution obscures behavior that a more realistic model 
would illuminate. 

In considering the current carried by cosmic rays 
through a cold ambient plasma parallel to a uniform mag- 
netic field, the calculations in B04 predict an instability 
in which the dominant mode, circularly polarized mag- 
netic fluctuations with wavenumber k\\ maxi grows at a 
rate 3w = -f max = v A k\\ ma x, ^ ~ 0, where 



k 



\\max 



eN C RV sh B\\ ^ LIqJcr 
~ 2Biio 



2N l m l v A 



(1) 



and va 



Bf ]0 /(i (N e m e + NiTUi) 



1/2 



is the Alfven ve- 



locity of the plasma. As mentioned above, this result 
is subject to the restriction that the fluctuation satisfy 
ijj <C S^i ; that is, that the growth rate of the instability 
must be much less than the nonrelativistic ion gyrofre- 
quency. 



2.2. Implementation 

Our simulations employ a modified version of the 
particle-in-cell code TRISTAN (jBunemanl [l99l . The 
code is updated to work in a high-performance comput- 
ing environment and u ses the charge-conse rving current 
deposition algorithm of lUmeda et al.l (|2003f) . We have in- 
corporated a fourth-order-accurate algorithm for updat- 
ing th e electromagnetic fields based on Maxwell's equa- 
tions ((Greenwood et al.ll2004D . and adapted the code to 
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TABLE 1 

Parameters and selected results of simulations 



Run 


Grid (A 2 ) 


Ni (A- 2 ) 


A se (A) 






ICR 


Bl/B|| 


T/ 'Ymax 


"~imax /^i 


A 


9000 x 2400 


25 


4 


0.4 


in 


50 


13.3 


0.96 


0.4 


A' 


3000 x 1600 


25 


4 


0.1 


in 


50 


17.6 


0.86 


0.4 


13 


4500 x 1800 


20 


6 


0,1 


hi 


50 


19.7 


0.89 


0.4 


C 


4000 x 1500 


25 


5 


0.3 


118 


10 


11.6 


0.78 


0.3 


c 


4000 x 1500 


25 


5 


0.3 


148 


10 


9.9 


0.93 


0.3 


D 


3200 x 1200 


25 


5 


0.4 


111 


10 


12.7 


0.86 


0.4 


E 


3200 x 1200 


25 


5 


0.4 


111 


25 


18.4 


0.90 


0.4 


F 


3200 x 1200 


25 


5 


0.4 


111 


200 


20.1 


0.91 


0.4 



Note. — Run A' reproduces the parameters of Run A on a spatial scale comparable to run B. Run C' differs from Run C only in that 
the cosmic-ray particles are not split into multiple simulated particles. JVj is the number of ions per cell (size A 2 ), A se = c/uj pe is the 
electron skin depth, v a h is the shock speed at which the cosmic-ray population drifts; ■ycR is t ne cosmic-ray Lorentz factor in the population 
rest frame. Ymra/Sii is the ratio of the predicted growth rate of the most unstable mode (with wavelength \ m ax) () to the nonrelativistic 
ion gyrofrequency. The peak measured perpendicular magnetic field and growth rate are given by B± and 7 in units of their respective 
reference quantities. Parameters common to all runs include an ion-to-electron mass ratio mi/m e = 50, a density ratio of ions to cosmic 
rays Ni/Nqr = 50, and an initial homogeneous magnetic field antiparallel to the cosmic-ray drift such that the Alfven speed is vj^ = 0.01c. 



a truly 2.5-dimensional simulation. 2 

The shortest relevant physical time- and length- 
scales are given by the electron plasma period u;~ = 

1/2. 

[m e eo/e 2 ./V e ] and electron skin depth \ se = c/oj pe , 
respectively. To ensure that physical small-scale elec- 
trostatic effects are distinct from those arising from the 
artificial discretization of space and time, for our simu- 
lations we choose parameters such that X se > 4A, where 
A is the length of one cell, and w^ 1 > 8.95t, where 5t 
is the timestep size. The wavelength X rnax = 2Tr/k\\ max 
of the predicted magnetic-field fluctuations is of order 
10 2 A se for all simulations (see Table [T]). We perform 
our simulations in a rectangular computational grid, with 
both the homogeneous magnetic field and the direction 
of cosmic-ray drift aligned with the longer dimension. 
Periodic boundary conditions in both directions allow 
particles and fields to "wrap around" to the far side of 
the grid, but this also imposes an upper limit to the size 
of structures that can be represented accurately. Having 
observed a migration of the dominant structures toward 
increasing size scales in N08, for one simulation (Run A) 
we chose a grid slightly larger than 20X max in the drift 
direction, and 5X max in the transverse direction. We ran 
additional simulations on smaller grids to explore the in- 
fluence of variations in our choice of parameters on the 
initial amplification of the magnetic field and the turbu- 
lent properties of the field and plasma. 

To minimize statistical noise associated with the par- 
ticles, we employ a combination of large per-cell particle 
counts and low-pass filtering of the electric currents aris- 
ing from their motion. We initialize the simulations with 
either 20 or 25 plasma ions per cell, and a density ratio of 
Ni/NcR = 50 plasma ions per cosmic-ray ion. However, 
in most runs we split each cosmic ray into ten simulated 
particles that preserve the charge-to-mass ratio but pro- 
vide better statistics, and the contribution of each parti- 
cle to currents and densities is weighted accordingly. 

Electrons, less massive than the ions by a factor of 
nii/m e = 50 in all tests, are physically one population 
but simulated as two populations. In a simulation for 
which Ni = 25, each cell contains 25 "full" electrons as 

2 As reported in N08, we have previously approximated two 
dimensions by simulating a volume with a thickness of only three 
grid cells. 
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Fig. 1. — Evolution of the magnetic field components in Run B 
displaying the growth and saturation of the instability predicted 
in B04. The magnetic field amplitude is in units of the initial 
homogeneous field, and time is in dimensionless units set by the 
calculated inverse growth rate 7^„ from that instability. 

well as five simulated particles, which add up to half of 
a split electron. The two simulated electron populations 
are initialized according to the same distribution func- 
tion, and behavioral comparisons of the full electrons and 
the split electrons in the simulations reveal no differences 
between the two populations, to within statistical limits. 

The initial temperature of the plasma is set artifically 
high, with u ej th — 0.01c, in order to mitigate Buneman- 
type electrostatic effects arising from the drift between 
the stationary population of ions and the slowly drifting 
electrons. Our simulations in N08 demonstrated that 
cooler plasmas were heated through such effects on a 
much shorter timescale than any turbulent magnetic- 
field amplification, but the resulting anisotropy of the 
ion distribution function in particular persisted on such 
timescales. A higher initial temperature leads to better 
preservation of isotropy against the intra-plasma drift. 
Additionally, using a density ratio Ni/Ncr of 50 instead 
of 3 (the value in N08) significantly reduces Vd- 

The cosmic rays are initialized to be mono-energetic 
and isotropic in the shock rest frame. This frame drifts in 
the — x-direction, antiparallel to the homogeneous mag- 
netic field, at a speed v s h = 0.4c or 0.3c. We set the 
magnetic field strength such that the Alfven speed va = 
0.01c, or equivalently that the plasma frequency exceeds 
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the Larmor frequency by a factor of uj pe /fl e « 14.14. 

2.3. Run-specific parameters 

In Table [1] we summarize the initial conditions and 
some observations from each of the simulations we dis- 
cuss. The default parameters are described below, fol- 
lowed by a list of each secondary run's purpose and any 
deviations from those parameters. 

Run A is our large-grid run for this paper, and is the 
basis for many subsequent "test" simulations. The cos- 
mic rays are isotropic and mono-energetic with Lorentz 
factor jcr = 50 in the CR rest frame, which drifts at 
v s h = 0.4c in the — ^-direction in the simulation frame. 
There are initially 25 plasma ions per cell, and each 
cosmic-ray ion is split into ten equal parts, for a net 
of five simulated particles per cell. For these parame- 
ters, Xmax — 444 A. In order to capture the non-linear 
stage when the dominant structure size exceeds this, we 
simulate on a grid of 9000A x 2400A, corresponding to 
20.26A maa; x 5.4A ma;r . 

Run A' uses the same physical parameters as Run A, 
but on a smaller grid, fitting 6.75X max X 3.6A max . 

Run A uses a large grid to ameliorate the limitations 
imposed by periodic boundary conditions, which is par- 
ticularly relevant during the later evolution when struc- 
tures in the magnetic field begin to grow quite large. 
However, for the purpose of analyzing the properties of 
the upstream medium during the growth and satura- 
tion, we simulate a smaller region with a greater detail 
level. Run B is our test with increased spatial resolution, 
in which the electron skin length has been increased to 
X se = 6A, and the particle density reduced to Ni = 20, 
while holding the other parameters fixed. The grid is 
G.73X max x 2.7A maa; . 

Runs C through F are tests with varying cosmic-ray 
Lorentz factors. All use the intermediate value for elec- 
tron skin length of A se = 5A. 

Run C has a reduced shock speed v s h = 0.3c. Its cos- 
mic rays also have a reduced Lorentz factor, jcr — 10- It 
is unique among our runs in that it does not split the cos- 
mic rays. The grid for Run C admits 5AX max x 2.03X max - 
Run C is identical to Run C but splits the cosmic rays. 

Runs D, E, and F use the usual shock speed v s h = 0.4c 
and cosmic-ray Lorentz factor jcr — 10, 25, and 200, 
respectively. Their grid measures 5.76X max x 2.16X max - 

3. RESULTS AND DISCUSSION 

The spatially-averaged amplitude of the turbulent 
components of the magnetic field, parallel and perpendic- 
ular to the initial homogeneous field Buq, are plotted as a 
function of time in Figure [1] for Run B, representing well 
the evolution in all runs. The unit of choice for time is the 
inverse predicted growth rate Jmax- When t^ rnax w 5, 
the field enters a period of exponential growth. As seen 
in Table [TJ the peak growth rate observed is often 10%- 
20% less than jmax', this small discrepancy may be pri- 
marily a thermal effect. As discussed in N08, the finite 
temperature of the plasma tends to reduce the maximum 
growth rate from its value in the cold-plasma limit. 

Initially, the growth of the magnetic field is stationary, 
in the form of non-propagating waves with wave-vectors 
oriented parallel to the cosmic-ray drift and similar in 
magnitude to kn max (Figure [2j The cosmic-ray Larmor 



radius rcR g 3> X max , so this is not the resonant inter- 
action of cosmic rays with self-generated Alfven waves, 
but is instead the non-resonant generation of an approx- 
imately purely growing wave mode. However, density 
fluctuations begin to appear while 5B <C B||o- The root- 
mean-square magnitude of the normalized density fluctu- 
ations (8Ni/Ni) during the linear growth stage is smaller 
than (SB/B^ a ), typically at a level of 20 to 50 % be- 
tween tj max = 5 and t^ rnax = 10. As the turbulent 
magnetic field becomes comparable in strength to the 
homogeneous field, the plasma begins to move and the 
cosmic rays slow down, reducing the current and satu- 
rating the magnetic-field amplification at a level of 10-20 

Any density and velocity fluctuations that arise may 
contribute to further deviation from a purely growing 
parallel mode, and may change the environment enough 
that further field amplification is not favored. We observe 
that when B± ~ _B|i , the growth rate is only gently de- 
creasing and the field still has approximately a parallel- 
mode structure, but when B± ~ 10-Biio, the spatial orga- 
nization of the field no longer resembles a parallel mode. 
As seen in Figure El the dominant length scale of the 
fluctuations grows and significantly exceeds X max . Sub- 
stantial variation appears along the direction transverse 
to the drift, reaching the maximum length scale allowed 
by the periodic boundary conditions of even our largest 
computational grid. The field growth comes to a stop, 
saturating generally between 10 and 20 times the ampli- 
tude of the initial homogeneous field. 

The reduction of the field growth rate accompanies the 
onset of turbulent plasma motion. Parcels of ambient 
plasma whose size is comparable to X max move and col- 
lide, giving rise to considerable density fluctuations. Fig- 
ure |4] (left) indicates that the plasma parcels can attain 
a considerable range of speeds in the directions trans- 
verse to the cosmic-ray drift. The strong perpendicular 
field and turbulent plasma allow bulk momentum trans- 
fer from the cosmic rays to the ambient plasma, so su- 
perimposed on the turbulent motions, the plasma begins 
to drift more and more rapidly in the direction of the 
cosmic rays, while the cosmic rays simultaneously slow 
down (Figure |H right), consistent with the quasi-linear 
predictions of lWinske and Lerovl (|1984f ) for non-resonant 
modes. In our simulations, the speeds converge such 
that the relative drift is roug hly an order of magnitude 
smalle r than its initial value. iRiquelme and Spitkovskvl 
(2009) suggest that convergence to a relative velocity ap- 
proximating the Alfven speed is sufficient to inhibit fur- 
ther magnetic-field amplification and thus that satura- 
tion may take place. Although defining the Alfven speed 
for a turbulent magnetic field is complicated, we find 
that by considering only the projection of the magnetic 
field onto the drift direction along with the local plasma 
density averaged over regions of a few X 2 se in area, we do 
observe a "mean" Alfven speed parallel to the cosmic- 
ray drift that differs from the post-saturation velocity 
separation only by a factor of order unity. 

In addition to slowing down in bulk, the instanta- 
neous rest-frame distribution of cosmic rays is modified, 
as shown in Figure [5] Some anisotropy is introduced, 
which may have consequences for the modeling both of 
radiation and particle acceleration. The anisotropy ap- 
pears in all simulations in a qualitatively similar man- 
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Fig. 2. — The perpendicular magnetic field in units of B^q (top) and density of plasma ions (bottom) at t-y m ax = 7.0 for Run A, when 
(5B/Bi|o) ~ 0.50. The parallel-mode structure is clearly visible in the magnetic field. Localized density fluctuations about the mean (25 
ions per cell) are smaller, with an amplitude of (SNi/Nj) xs 0.08. 
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Fig. 3. — As in Figure [2] the perpendicular magnetic field and density of plasma at t"f max = 14.0 in Run A. The plasma density 
fluctuations are strong and well correlated with the magnetic field structures during the non-linear stage of magnetic field growth. The 
structures move rapidly on oblique trajectories at this time in the simulation (see Figure [4}. 



ner, though the rate at which it develops (as well as the 
extent to which the CR bulk speed changes) decreases 
with increasing jcr- The fastest development of rest- 
frame anisotropy occurs for cosmic rays whose trajec- 
tory is transverse to the drift direction in the simula- 
tion. During the linear stage of magnetic-field growth, 
the stationary parallel- wave mode with wavelength X max 
will be favorable for resonant scattering of cosmic rays 
whose pitch angle /i is given by Xmax^CR/^ = v \\ = 

1 /2 

fxc (l — 7^fl) • For the parameters i n our simulations 
and i n agreement with the work of ILuo and Melrose] 
/x is of order 10 2 , corresponding to those cos- 
mic rays in which the aforementioned anisotropy is first 
seen. Additional scattering and acceleration of cosmic 



rays may occur in the electric fields associated with mag- 
netic fields being transported in the turbulent flow dur- 
ing the nonlinear stage. Second-order Fermi acceleration 
is possible in such a turbulent environment. The elec- 
tric fields we observe during saturation are typically an 
order of magnitude smaller than the magnetic fields. A 
Lorentz transformation into the global plasma rest frame 
reduces the electric fields only from about t^ max = 15 
onward, in agreement with the behavior of the plasma 
bulk speed. If we transform the EM fields into the local 
plasma rest frame, thereby accounting for its turbulent 
motion, we see a reduction of the electric fields beginning 
at t^max ~ 10. After about t^ max = 15, less than 20% 
of the original electric field remains, demonstrating that 
the dominant fraction of the electric field in our Simula- 
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Fig. 4. — Left: From Run B, the transverse in-plane velocity of the plasma in regions of size (0.45A ma;E ) 2 . Two regions have been selected 
and superimposed in lighter colors for clarity. Right: In the same run, momentum transfer from the cosmic rays to the background plasma 
removes their relative drift, limiting the amplification of the magnetic field. The "overshoot" at t^ ma x ~ 16 appears in the small-grid runs 
but is not present in Run A, with the larger grid, so it is likely a consequence of the periodic boundary conditions 



tion arises from the transport of magnetic fields. This 
indicates that pure elastic scattering may be a poor as- 
sumption, if one considers the diffusive transport of high- 
energy cosmic rays in the upstream region of young SNR. 
Runs C-F with varying cosmic-ray Lorentz factors illus- 
trated principally that "stiffer" cosmic rays took slightly 
longer to begin slowing down and were marginally more 
resistant to the onset of rest-frame anisotropy, but there 
was no qualitative difference in the linear growth of the 
magnetic field nor in the saturation mechanism. 

The interactions between the shock itself and any 
spatial variations not only in plasma density but also 
in velocity are known to introduce vorticity in the 
downstream region that may amplify the magn etic field 
(jGiacalone and Jokipiil 120071: llnoue et all 12009ft . They 
may also have a non-negligible impact on the injec- 
tion process or the scattering properties of the down- 
stream medium, and consequ ently the resulting distri- 
bution of energetic particles (jMalkov and Drurvl 120011) . 
If the injection mechanism depends primarily on highly 
localized processes, as opposed to being a product of 
the average conditions over large spatiotemporal scales, 
then the presence of appropriately distributed inhomo- 
geneities might reduce the availability of readily accel- 
erated particles in one region, with no guarantee of 
a corresponding increase in a neighboring region. In 
particular, if injection is ultimately more efficient at a 
parallel shock than a perpendicular shock, the devel- 
opment of a turbulent magnetic field in the formerly 
parallel region via this streaming instability would re- 
duce the overall quasi-parallel fraction of the shock's 
surface if comparable turbulence is not established in 
the initially quasi-perpendicular regions. Beyond injec- 
tion, the cosmic-ray transport properties upstream and 
downstream contribute significantly to the acceleration 
efficie ncy and the maximum energy cosmic rays can at- 
tain. IMarcowith et alj (|2006f ) investigated the sensitivity 
of the cosmic-ray diffusion coefficient to the character of 
the turbulence, and found that, in particular, anisotropic 
turbulence may lead to somewhat higher maximum par- 
ticle energy. 

In addition to possible implications for the spectrum of 
accelerated particles, the character of the upstream mag- 
netic turbulence also eventually contributes to the emis- 
sion properties of the downstream medium. The pres- 



ence of anisotropy in the turbulence after being shock- 
compressed may lead to a larger degree of polarization 
in the radio synchrotron emission tha n might otherwise 
be expected (jStroman and Pohl 200^). 

The turbulent motion in the upstream medium may 
be significant for self-consistent models of cosmic-ray ac- 
celeration. In our simulations, plasma parcels attain ve- 
locities that exceed the initial Alfven speed and collide, 
producing "fronts" whose effects may include further lo- 
calized field amplification or, perhaps quite importantly, 
compression and heating of the upstream medium. We 
observe that heating of the ambient plasma occurs cospa- 
tially with the interfaces between rapidly moving plasma 
parcels and the slower regions into which they are press- 
ing. In Run B, for instance, the average kinetic en- 



ergy per ion is of order 10 m^c at tj n 



3.8, be- 



fore the linear growth has begun, and is slightly higher 
but still of that order at t~f max — 7.6, during linear 
growth. But at t^ max = 11.5, the beginning of satu- 
ration, the mean kinetic energy per ion has increased to 
~ 5 x 10 _4 miC 2 , and at t^ max = 15.3, it has climbed to 
10~ 2 miC 2 . The common and convenient assumption that 
the entire amplification and saturation process proceeds 
in a low-temperature regime must be carefully consid- 
ered if conditions allow for turbulence to develo p rapidly 
in the precursor, but as lVladimirov et ail (|2008ft observe, 
the substantial heating expected from the dissipation of 
even a small fraction of upstream turbulence may signif- 
icantly increase the rate at which particles are injected 
into the acceleration process. 

The reduction of the cosmic-ray current (by the con- 
verging drift velocities of their population and the am- 
bient plasma) serves to underscore the need to treat 
these particles kinetically. The back-reaction from the 
turbulent plasma and amplified magnetic field is non- 
negligible, with the most prominent consequence being 
saturation of the field growth at a lower amplitude than 
might be expected from a constant cosmic-ray current 
(jOhira et all (2009ft . 

Finally, it is worth noting that although the initial 
dominant mode of magnetic-field amplification agrees 
with the calculations of B04, the later evolution and sat- 
uration mechanism have much in common with the sim- 
ulations performed in N08, in which the choice of simu- 
lation parameters exposed an oblique filamentary mode 
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as the dominant instability. It seems it doesn't matter 
what form the initial linear instability takes. 

An important consideration in interpreting our find- 
ings is that the ability of our simulations to accurately 
portray the evolution of the magnetic field and the 
plasma properties at very late times, beyond the satura- 
tion, is limited. In addition to the limitation imposed by 
the finite size of the simulation domain, the assumptions 
of environmental homogeneity become increasingly un- 
sound. Changes in the plasma flow speed i mply density 
adjus tments in order to conserve mass flux (jOhira et alJ 
2009). Also, as the shock approaches, spatial gradients in 
the cosmic-ray precursor and ambient plasma properties 
play an increasingly important role. Furthermore, the 
time in which upstream amplification can occur before 
the system is overtaken by the shock is only ~ ^Imax 
for our parameters, assuming efficient Bohm-type diffu- 
sion so that the shock sweep- up timescale is crcR a l'&v 1 sh 
(where rcR g is the cosmic-ray gyroradius) (N08). Thus, 
field growth is inevitably limited even in the absence of 
saturation. 

It bears mentioning that to a first approximation, the 
effects of a changing environment will be limited. As 
mentioned in §2.11 the parameters of the streaming insta- 
bility are agnostic to the details of the cosmic-ray distri- 
bution. Thus, although an approaching precursor implies 
a continual softening of the local cosmic-ray energy spec- 
trum due to the arrival of particles over an increasingly 
inclusive energy range, it is only the increased overall 
flux that matters. Until surplus high-energy electrons 
are included, the net cosmic-ray current will increase, 
whose effect according to Equation [1] will be to amplify 
the parallel mode described therein on shorter time- and 
length-scales, potentially well separated from the large 
scales attained by the turbulence generated earlier. Even 
this effect is difficult to predict with confidence: not only 
will cosmic rays that arrive during saturation encounter 
a turbulent upstream medium, but may themselves have 
been influenced by prior unstable evolution nearer to the 
shock. In any case, since the precursor length scale is 
proportional to the diffusion coefficient, and hence to the 
energy for Bohm diffusion, cosmic ray particles of signifi- 
cantly lower energy (and higher flux) are brought in only 
shortly prior to capture by the shock. 

Although the non-resonant streaming instability 
is concerned only with the current t hey carry 
(and, as demon s trated by iPelletier et al.l (|2006f ) and 
ILuo and Melrose! ((2009), is largely produced by the re- 
turn current in the plasma), this or other instabilities 
thus encountered may introduce anisotropy and inhomo- 
geneity into the local cosmic-ray distribution, with po- 
tentially non-trivial consequences. As increased compu- 
tational capacity and improved models become available, 
we can expect simulations representing a significant ex- 
tent of the precursor to yield valuable insight into these 
interrelated and nonlinear processes occurring there. 

4. SUMMARY AND CONCLUSIONS 

We have simulated the turbulent amplification of the 
interstellar magnetic field upstream of a nonrelativis- 
tic shock using a kinetic 2.5-dimensional particle- in-cell 
code. 

We observe: 
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Fig. 5. — The phase-space distribution of cosmic rays in the in- 
stantaneous rest frame of the population at time t-fmax ~ 15.7 in 
Run B, compared with the initial distribution (thin semicircle). In 
addition to revealing the anisotropy in the forms of stretching and 
broadening that began during the linear magnetic-field growth, the 
instantaneous rest frame differs from the initial frame by approxi- 
mately 0.18c. 

• That the non-resonant streaming instability (B04) 
is seen initially for our choice of parameters, but 
with fluctuations appearing in the plasma density; 

• That the plasma quickly evolves and saps the bulk 
momentum from the drifting cosmic rays; 

• The reduction of relative drift effectively removes 
the cosmic-ray current and saturates the magnetic- 
field amplification to ~ 20-Bi|o; 

• This saturation mechanism is independent of the 
initial linear instability, occurring as a result of the 
back-reaction on the cosmic rays; 

• Strong transverse plasma motions arise in conjunc- 
tion with the turbulent magnetic field exceeding 
the homogeneous background field; 

• Interactions between differently-moving plasma 
parcels give rise to significant fluctuations in den- 
sity and temperature of the upstream interstellar 
medium; 

• The plasma begins to drift, suggestive of a cosmic 
ray-modified shock, until its speed approximately 
matches that of the cosmic-ray population. 

• The phase-space distribution of cosmic rays un- 
dergoes significant changes: the development of 
anisotropy within the population rest frame begin- 
ning during the linear growth of the magnetic field, 
and decrease in bulk speed associated temporally 
with the saturation of the magnetic field; 

• Evolution past this point is beyond the reach of our 
current simulation method, as it requires additional 
input to reflect the changing external environment. 

The growth of non-resonant streaming instabilities in 
the cosmic-ray precursor of supernova shocks may con- 
tribute to significant turbulent amplification of the mag- 
netic field upstream of the shock. However, the ampli- 
tudes achieved via any one process are likely to be some- 
what less than suggested in the artificial case of constant 
cosmic-ray current. Simulations that incorporate incre- 
mentally more of the environmental changes expected in 
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a real precursor scenario can be expected as computa- 
tional capacity increases, although it may take innova- 
tions beyond the foreseeable future to make a realistic 
particle-in-cell app roach tractable for such a simulation 
(jVladimirov et al.ll2008h . 
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